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This  report  summarizes  the  efforts  at  UT  Dallas  re  evaluate  the  .  re 
measuring  capability  of  the  AFGL  accelerometers  (MKSA'si  mounted  on  the 
Atmosphere  Explorer  satellites.  The  so'-  two  re  available  let  thi:.  purpose 
AFGL  had  some  I  imitations  that  prevented  its  use  lor  our  iuirpo-.es ,  am!  mm-h 
contract  time  was  devoted  to  producing  -isnu  il<  I ■  e  i  i  w.i  r-  pa.  I  ij-e  .  m; 
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Both  AE-D  and  AE-E  have  another  device  i NATE)  on  hoard  that  measure 
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same  wind  component  as  MESA.  Our  limited  comparisons  of  these  measurements 
show  that  they  are  in  reasonable  agreement.  We  suggest  that  a  substantial 
amount  of  MESA  wind  data  from  these  two  satellites  be  processed  in  order  to 
make  a  careful  comparison  with  the  NATE. 

On  AE-C  NATE  does  not  measure  the  horizontal  wind  component,  and  we 
suggest  that  the  MESA  wind  data  for  this  satellite  be  Drocessed  in  its 
entirety  (for  the  elliptic  orbit  phase)  and  be  placed  in  the  UA  file. 
Fortunately,  the  AE-C  MESA  data  are  the  cleanest  of  all. 
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INTRODUCTION 


Tne  AFGL  Miniature  Electrostatic  Accelerometer  (MESA)  on  the  Atmosphere 
Explorer  (AE)  satellites  was  designed  to  provide  neutral  atmospheric  density 
measurements  by  determining  satellite  deceleration  caused  by  ae  odynamic 
drag.  Each  spacecraft  contained  three  single-axis  sensors.  The  first  is 
aligned  with  the  Z  (spin)  axis,  while  the  other  two  sensors  (designated  XT' 
and  YX)  are  located  in  the  X-Y  plane  at  a  45°  angle  on  either  side  of  the 
+X  axis.  The  sensors  have  the  capability  of  being  commanded  into  any  of 
three  sensitivity  ranges,  with  the  least  sensitive  being  only  used  for 
satellite  thrust  monitoring.  The  focus  of  our  effort  has  been  the  utiliza¬ 
tion  of  the  data  from  these  three  sensors  to  determine  the  velocitv  of  the 
component  of  neutral  atmospheric  winds  lying  along  the  spacecraft  Z  axis 
during  those  periods  when  spacecraft  altitude  was  low  enough  ('  L80  km)  for 
a: curate  measurements. 

Thu  first  section  of  the  report  discusses  the  software  package  that  was 
put  together  at  UTD  to  calculate  the  winds  from  the  MESA  telemetry  data. 

This  is  followed  by  a  discussion  ol  the  application  of  this  software  to 
actual  MESA  data.  Results  are  then  plotted  and  compared  with  NATE  data,  and 
the  outputs  of  ditterent  filters  are  examined.  Finally,  some  recommendations 
tor  MESA  data  analysis  are  presented. 


SOFTWARE 


Hie  first  step  was  the  examination  ol  the  results  of  the  AFGL  MESA  Data 
Reduction  System  (DIM),  which  im  l tided  a  preliminary  wind  estimate.  It  was 
initial Lv  apparent  that  the  coat r ibut ion  to  tin*  winds  of  atmospheric  corota— 


tion  >00  m/see )  was  not  included.  In 
suddon  Large  discontinuities  (.see  Figun 
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addition,  the  computed  values  showed 
!)•  The  cause  of  these  effects 
ruble  to  investigate  the  effectiveness 
'tom  the  d.ita.  For  these  reasons 
program  to  ui  ilv.te  the  MFSA  data  and 
pack  ipe  deve,oi>ed  is  described  in 
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To  simplify  calculations  and  to  achieve  maximum  accuracy  ii  was  decided 
to  use  only  data  t rom  despun  passes.  In  addition,  the  program  cues  not  re¬ 
quire  any  input  from  the  MESA  data  reduction  file  system  developed  at  the 
AFGL,  using  only  the  TM  data  base  and  orbital  and  attitude  information  from 
the  AE  Central  Computer  Facility  at  Goddard  Space  Flight  Center.  The  pro¬ 
gram  converts  the  TM  data  to  MESA  raw  data,  applies  a  temperature  correction 
to  the  data,  utilizes  a  filter  to  eliminate  noise  due  to  the  satellite's 
momentum  wheel  assembLy,  and  finally  calculates  the  velocity  of  the  neutral 
wind  along  the  Z-axis  of  the  satellite,  producing  a  line  printer  plot  of  Z 
axis  wind  versus  G.M.T.  The  Z  axis  wind  was  calculated  from  the  formula: 
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where  W„  is  the  Z  axis  wind,  a„  is  the  Z  axis  acceleration,  V  '  is  the  satel- 
..  Z  s 

lite  velocity  in  the  earth's  rotating  frame  of  reference,  a+  is  the  total 
satellite  acceleration,  and  a  is  the  angle  between  the  Z  axis  and  the  satellite 
velocity  vector.  Total  satellite  acceleration  may  be  determined  either  by  using 
the  output  from  aLL  three  sensors  or  bv  utilizing  only  the  XY  and/or  YX  sensor 
and  correcting  for  the  angle  between  the  sensor  and  the  spacecraft  velocity 
vector. 


Control  and  Input  Data 

basic  program  control  is  simple.  The  major  required  parameters  are 
satellite  idem  i  t  i  cat  ion ,  start  and  stop  dates  and  times,  sensor(s)  to  be 
used  to  compute  total  drag,  reference  frame  for  wind  calculation  fcorotating 
or  rion-eorot.it  tug)  ,  and  number  of  data  points  used  in  the  filter.  The 
program  then  proceeds  to  extract  from  the  AE  data  base  all  data  necessarv, 
including  raw  I'M  data,  sensor  temperature,  range  indication,  overflow  flags, 
and  osition,  attitude  arid  velocity  information. 

I  omper.iture  Correction  and  Filterin', 

'  c.  !.:  data  is  I  irst  temperature  corrected  and  then  led  into  a  buffer 
Iron,  where  it  is  :  non  passed  through  a  i  i Iter  to  remove  momentum  wheel  and 
nut  i!  a  no  iso.  one  .«!  the  I  liters  used 
TI'''  leg  use  .in  'lie  Midi  A  ivstem. 

,;Sei 1  \e iiorid  i  x  . 


is  the  "l.-lo"  filter  developed  at 


Wind  Determination 


Data  from  above  280  km  altitude  is  used  to  establish  a  xero  accelera- 
tion  baseline  for  sensor  response.  Then  using  the  XY  sensor, 
the  YX  sensor,  or  all  three  sensors,  the  total  drag  induced  deceleration  of 
the  vehicle  is  determined.  Using  the  known  attitude  and  velocity  01  the 
satellite.  Equation  1  is  then  applied  to  yieLd  the  7.  axis  wind. 

Overranging 

When  one  of  the  sensors  is  in  its  most  sensitive  range  it  can  uappen 

— 

that  the  acceleration  along  that  axis  may  exceed  the  maximum  (2  x  10  g) 
for  that  range.  When  this  occurs  an  overrange  flag  is  signaled  in  the 
telemetry  for  that  sensor.  When  the  XY  or  YX  sensor  is  being  used  to  deter¬ 
mine  total  drag,  an  overrange  condition  causes  that  data  point  to  be  dropped 
In  intervals  where  the  Z  axis  sensor  overranges  no  useful  wind  information 
can  be  obtained  from  the  present  software,  but  it  is  likelv  that  data  could 
be  recovered  in  these  regions  with  a  little  more  effort.  This  ususallv 
occurs  only  below  140  km  altitude. 

Output 

The  program  output  consists  of  four  data  files.  The  1  irst  consists  of 
processed  analog  data.  The  second  contains  processed  digital  data  (this  is 
the  main  output  file).  The  third  contains  a  Line  printer  plot  of  the  7.  axis 
wind.  The  Last  is  a  file  indicating  data  dropped  due  t  >  errors  or  over- 
ranging. 

PROGRAM  APPLICATIONS 

Since  there  is  no  existing  wind  data  base  from  NATE  lor  AK-G  this  salel 
Lite  performance  is  most  critical  from  an  overall  AE-se ience  eons  iderat i on . 
In  the  course  of  data  reduction  several  constraints  upon  the  usage  of  the 
Winds  program  in  regard  to  AE-C  data  become  apparent. 

One  problem  exists  regarding  the  calibration  of  the  XY  and  YX  sensors. 
If  during  the  course  of  a  low  altitude  pass  one  uses  the  aiceler.it  on  com¬ 
puted  from  one  of  these  sensors  and  the  orientation  angle  O  that  onset 
relative  to  the  direction  ot  satellite  velocity  one  '.in  obtain  a  value  tor 
the  total  drag  induced  aci  ••  1  er.it  ion  ot  the  v  •  ■  1 1  i  •  l « • .  It  i  he.  is  done  toi 


t>uch  of  tile  two  spin  plane  sensors  one  consistently  finds  that  they  are  in 
disagreement,  with  the  XY  sensor  giving  values  approximately  10  to  20% 
larger  than  those  from  the  YX  sensor.  This  seems  to  indicate  the  need  for 
an  in-! light  calibration,  and  this  can  perhaps  be  achieved  bv  til iliaing  the 
density  measurements  from  the  onboard  neutral  mass  spectrometers. 

A  second  concern  is  that  of  the  sensitivity  ranges  of  the  sensors.  The 
standard  mode  of  operation  appears  to  have  been  to  leave  the  YX  and  Z 
sensors  in  their  most  sensitive  range  (range  3)  while  leaving  tae  XY  sensor 
in  range  2.  While  this  practice  yields  the  maximum  amount  of  information  in 
regard  to  total  atmospheric  density,  it  requires  that  the  XY  and  YX  sensors 
be  intercal Lbrated  to  produce  a  smooth  wind  profile,  since  each  is  used  for 
part  of  each  pass.  lhe  YX  sensor  saturates  and  overranges  below  150  km 
altitude,  forcing  us  to  use  only  the  XY  sensor  for  total  drag  determination 
below  this  altitude. 

In  the  case  ot  the  AE-P  and  AE-E  satellites  it  was  found  that  the  XY 
and  YX  were  also  left  on  ranges  2  and  3  respectively.  The  Z  sensor,  however, 
was  not  set  on  Lhe  most  sensitive  range  but  instead  on  range  2.  This  de¬ 
grades  the  accuracy  of  wind  determination  using  these  two  spacecraft,  but  not 
by  a  large  t actor,  because  the  signal  to  noise  level  is  fairly  low  so  that 
the  extra  gain  (x  20)  going  from  range  2  to  range  3  does  not  help  too  much 
(i.e.  the  s i gna 1  is  large  enough  to  be  seen  on  range  2.) 


WIND  RESULTS 


Hecau.se  ol  the  short  t  iine  available  to  examine  MESA  wind  data  after  the 
new  otiware  was  checked  out  the  conclusions  reached  here  are  somewhat  tenta¬ 
tive.  A  total  oi  approximately  40  orbits  from  all  three  satellites  provide 
the  basis  for  this  report. 

lhe  Ah-D  data  plotted  in  Figure  L  show  a  comparison  of  the  east-west 
wind  component  as  determined  i-v  tlu-  program  developed  here,  by  the  AFGL  pro¬ 
gram,  and  bv  the  NATE  instrument.  Except  for  the  regions  where  the  AFGL 
result  has  an  apparent  sign  difficulty,  the  three  sets  of  data  are  in  general 
agreement.  i'he  smoother  appearance  of  the  AFGL  output  compared  to  the  LTD 
point  ;  arises  because  ol  much  heavier  filtering  (a  two-minute  filter  compared 
lo  a  tu  second  tiller).  Some  oi  tlu*  A E-P  MESA  data  is  quieter  tiian  that 
shewn  lure,  and  s.sm*  is  noisier,  but  there  is  much  useful  data  even  though  the 
•  .ixis  , »cri‘  1  t'ronu' !  •  •  r  *\ a:,  on  tlu*  miiliilo  stMis  i  t  i  v  i  iv  r.uiii.r. 
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Perigee  was  near  82°  north  latitude  an  orbit  D-437  so  that  the  wind 
pattern  may  be  reasonable  in  spite  of  its  structured  appearance. 

A  more  detailed  comparison  of  the  various  filter  outputs  is  shown  in 
Figure  2,  which  also  uses  data  for  orbit  D— *57  (note  the  change  in  sign  of 
the  wind  coordinate).  Data  points  are  plotted  everv  10  seconds.  All  the 
curves  have  pretty  much  the  same  general  characteristics  except  for  the 
slightly  higher  noise  level  with  the  narrower  filters.  The  simple  running 
filter  seems  to  be  as  effective  as  the  much  more  sophisticated  digital 
filter. 

There  is  no  wind  data  with  whicli  to  compare  on  AE-C,  but  the  results 
obtained  from  the  sensors  using  current  software  are  consistent  with  those 
expected  on  a  theoretical  basis.  Figure  3  presents  two  typical  equatorial 
passes  recorded  nineteen  days  apart  on  the  AE-C  satellite.  Note  that  in 
both  cases  there  exist  data  gaps  near  perigee  due  to  the  saturation  of  the 
2  sensor  (not  all  points  are  saturated,  only  those  with  large  noise  spikes). 

In  addition,  the  effects  of  atmospheric  corotatlon  have  been  removed,  coro¬ 
tation  introducing  an  easterly  bias  of  approximately  450  m/sec  near  the 
equator.  The  perigees  of  orbits  2830  and  3086  occur  at  a  Solar  Local  Time 
(SLT)  of  about  16:40  and  12:00  respectively.  Models  of  the  wind  velocity 
determined  from  atmospheric  pressure  gradients  predict,  zonal  wind  velocities 
of  about  100-150  m/sec  in  these  regions,  with  the  wind  at  L 2:00  SL'!  being 
westward  and  the  wind  at  16:00  SLT  being  eastward.  Indeed,  the  data  do 
reveal  a  westward  to  eastward  shift  in  the  wind  velocity  near  perigee  with 
the  magnitude  of  the  change  being  about  110  m/sec.  The  noise  level  of  the 
rest  of  the  AE-C  data  examined  is  of  comparable  magnitude.  In  general, 
this  noise  Level  is  probably  smaller  than  the  absolute  uncertainty  in  the 
winds  arising  from  errors  in  the  spacecraft  attitude.  The  results  eertainlv 
appear  to  be  of  useful  quality  for  AE-science  purposes. 

A  comparison  of  the  NATE  wind  data  from  AE-E  with  that  from  Ml  SA  is  shown 
in  Figure  4  for  orbit  385.  There  is  a  general  consistency  between  the  two 
results  though  there  is  an  offset  of  40  to  30  m/s  near  perigee.  This  dis¬ 
crepancy  is  not  a  great  deal  greater  than  the  uncertainty  in  the  mounting 
position  of  the  devices  within  the  spacecraft.  From  a  comparison  of  u 
large  data  base  of  NATE  winds  with  the  ion  drift  meter  results  from  AE-E,  we 
also  find  an  offset  between  the  devices  of  the  order  of  '.<>  m/s.  This  offset 
is  slowly  varying  in  time  (weeks  to  months)  and  is  noi  vet  understood. 
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SUMMARY  AND  RECOMMENDATIONS 


In  summary,  the  technique  of  using  accelerometers  <>n  a  low  altitude 
satellite  to  determine  the  transverse  neutral  wind  appe.tts  to  provide  v.il  id 
measurements.  Though  a  relatively  high  noise  level  is  pesenl  on  the  AH 
satellites,  the  current  MESA  sensors  are  adequate  to  provide  a  reasonable 
wind  determination,  particularly  when  the  Z  axis  sensor  is  in  its  most 
sensitive  range,  as  in  the  case  ot  AE-C.  Al.-I)  and  All -K  orbit. t,  wiiieh  have 
a  slightly  higher  noise  level,  are  nevertheless  valuable  in  that  NATE  wind 
data  provides  a  comparison  between  different  techniques  and  ran  he  used  to 

build  confidence  in  the  AE-C  orbits.  This  data  should  be  completely 
analyzed  and  the  results  placed  in  the  Unified  Abstract  (UA)  files  for 
reference  by  all  interested  researchers.  A  word  assignment  lot  MESA 
wind  already  exists  in  the  UA  files. 

We  did  not  expend  any  effort  to  examine  the  usefulness  of  the  spinning 
data,  but  it  is  probably  worth  some  effort  to  establish  whether  this  data 
is  of  any  value  also. 

It  may  be  that  the  MESA  data  from  both  AE-D  and  AE-K  should  also  be 
completely  analyzed  and  placed  in  the  UA  files.  Certainly  a  relatively 
thorough  comparison  with  the  NATE  data  should  be  made,  if  for  no  other 
reason  than  to  determine  the  proper  confidence  to  he  placed  in  the  technique 
for  the  AE-C  wind  users.  This  is  probably  a  matter  that  might  be  discussed 
with  the  AE-science  team. 

At  UTD  we  have  developed  a  data  analysis  program  monitor  in  connection 
with  the  processing  of  the  RPA  and  MIMS  data  from  All.  il  this  scheme  were 
implemented  for  MESA,  the  further  effort  involved  to  complete  the  data 
analysis,  place  the  data  in  the  UA  file,  and  provide  a  lilm  plot  ot  each 
orbit  processed  would  be  minimal.  The  program  would  have  to  run  in  the 
MESA  Account.  If  this  plan  is  to  be  given  consideration,  it  should  he 
implemented  soon,  because  of  the  impending  loss  ol  the  h  computer  for  AE 
data  process -ng. 
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APPENDIX 


Neutral  Wind  Measurement 

We  desire  to  use  the  measured  drag  force  K  on  a  satellite  to  <K  t ermine 
the  neutral  gas  winds  (W) .  We  can  write  (in  the  satellite  t  rame) : 

V  =  -V  +  <o  x  R  +  W  (  I  i 

n  s 

where  V  =  total  neutral  gas  velocity  in  the  satellite  I rame 
n 

Vg  =  satellite  velocity 

to  =  earth's  angular  velocity 

R  =  satellite  position  vector 

W  =  neutral  wind  in  earth's  rotating  frame  of  reference. 


The  drag  force  on  the  satellite  is: 

F  =  1/2  p  ACd  |Vn|Vn  (2. 

where  p  =  atmospheric  density 

0^=  drag  coefficient  in  direction  of 
A  =  presentation  area  normal  to  ot  satellite. 

Assume  the  satellite  is  oriented  with  its  X  axis  approximately  in  the  ram 
direction  and  with  its  Z  axis  horizontal.  T ten  the  acceleration  components  are: 


a  =  F  •  Z/M  (3) 

7. 

=  7  •  X/M  (.) 


where  M  =  satellite  mass. 
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We  also  note  that  V  is  parallel  to  a  where  a,  =  a  X  +  a  Y  +  a  Z. 

n  +  +  x  y  2 

Making  the  substitution  V  '  =  V  -  U)  x  R  we  can  rearrange  (1): 

W  =  V  +  V  '  (5) 

n  s 

or  looking  at  just  the  7.  component 


W  =  V  +  V  • 

7.  nz  s  z 


Using  e‘(.  (1)  and  (4)  we  can  show 


(6) 


a 

z 


a 

x 


V 

=  nz 

V 

nx 


(7) 


l.'s  i  iu\  th>*  def  ini  t  ion  < > !  V  : 

nx 


V  -  V  •  X 
nx  ti 


X 


i 

i 


v 

n  ' 


(8> 


l  .  im:  ei; .  ( 
express  ion 


2)  i  •>  el  i in i  nate 

for  V  into  eq 
nx 


V  in  eq.  (8)  we 
'  n 

(  7 )  and  ge  t : 


then  substitute  this 


V 

n.‘ 


:  M 


1/2 


PA  C 


1) 


i  9) 


Thun  the  expression  lor  tlie 


component  of  the  neutral  wind  becomes: 


i  ( 


M 


.1/7 

•)  +  V 


(  10) 


'  depends  only  on  satellite  position,  velocity  and  the  earth's 
rotation.  These  parameters  can  be  obtained  from  the  orbital-attitude  data  of  the 
A.E.  computer  facility.  However,  an  independent  value  of  p  is  also  roqu i red . 

Alternatively,  using  equations  (6),  (7),  and  (8)  we  can  eliminate  the 
dependence  of  on  the  magnitude  of  p  or  C^,  and  obtain: 


W  =  V 
z  1  n  1 


Thus  the  fact  that  and  A  depend  on  the  angle  of  attack  is  cancelled  out 
if  Equation  (11)  is  used.  But  even  if  Equation  (10)  is  used,  it  is  not  likely 
that  CD  will  change  very  greatly  over  the  small  angles  of  attack  normally 
encountered  (<10°). 


Making  the  approximation  |V  j  =  |v  'I  Equation  (11)  becomes: 


n  ■  s 


W  =  |V  +  V  '  .  z 

3  la+| 


-  <a,  +  a.  cosa> 

W.  =  I  V  •  |  — ± ± - 3  W/ 

'■1  I 

'  +  ' 


where  a  «  angle  between  V  '  and  Z. 


Equation  (13)  is  the  same  as  Equation  (!)  in  the  text,  and  is  used  in  the 

UTD  winds  program  to  determine  W^.  Without  an  accurate  independent  measurement 

of  p  no  determination  can  be  made  of  the  other  components  of  W. 

Estimate  of  Error.  The  validity  of  the  approximation  Iv  1  =  1 V  * f  tnav  be 

1  n  1  s 

examined  under  typical  conditions. 

Using  (11)  and  (13) 


W  '  +  ( 1 V  J 

z  '  n 


we  know  ( ' V  !  -  ! V  ' i )  <  ' W ! 

n  s  — 


and  typical  lv  a  /'a  i  <  0.08 

z  + 


f 


So  for  a  200  m/s  wind  the  error  term  in  (14)  is  <  16  m/s. 

A  secondary  source  of  error  is  introduced  in  that  | a+ |  is  typically 
determined  by  using  only  one  sensor  and  making  the  assumption  that  the 
neutral  wind  velocity  is  anti-parallel  to  the  satellite  velocity,  thus 
we  can  write: 


a  /  (X 
x 


ax  >  <* 


-  V 

s 


|vs’| 


(15) 


For  | W |  =  200  m/sec  this  leads  to  a  maximum  2.17,  error  in  |a+|.  Using 
(12)  this  impliei  an  error  <  16  m/s.  Thus,  these  two  sources  of  error  are 
about  equal  in  magnitude. 


